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£T) ' Abstract 

o ' 

\Q . The 3d gluodynamics which governs the large T quark gluon plasma 

is studied in the framework of the field correlator method. Field cor- 
relators and spacial string tension are derived through the gluelump 
Green's functions. The glueball spectrum is calculated both in C = — 1 



Qh! as well as in C = +1 sectors, and multigluon bound states in the form 



of "gluon rings" and "gluon stars" are computed explicitly. Good 
overall agreement with available lattice data is observed. 



x. 

1 Introduction 



As was suggested fifteen years ago [TJ [2J , in QCD magnetic fields survive the 
temperature phase transition and keep their magnitude almost intact up to 
T = 1.5 T c , as was checked by lattice measurements of field correlators [3 . 

At larger T, T > 2T C , magnetic fields as measured by the spacial string 
tension a s (T) start to grow quadratically, cr s {T) ~ T 2 |HEj, which signifies 
the advent of a new regime, called the dimensional reduction. In this regime 
the temporal direction is squeezed, At = ~ — > 0, while higher Matsubara 
frequencies are suppressed, so the effective dynamics is reduced to 3d gluo- 
dynamics [Hj-jni- The latter is a confining theory and usually one relies on 
lattice 3d calculations to obtain physical quantities, such as cx s (T) |HIZ|. 
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Instead we suggest here to use the Field Correlator Method [TU] , (see [TT] 
for a review) to study the 3d nonperturbative dynamics and in particular to 
compute field correlators, string tension and lowest mass excitations in 3d 
gluodynamics analytically. Knowing cx s (T) one can compute other physical 
quantities of interest, e.g. in the 3d screening masses of mesons and 
of the lowest glueball were obtained in good agreement with existing lattice 
data. 

To find cr s (T) one needs to know the nonperturbative correlator of mag- 
netic fields D H (x) Pig. The latter can be studied analytically using the new 
method suggested in [T3], where it was expressed in terms of the gluelump 
Green's functions. In case of Ad QCD the gluelumps were studied in [T3] , and 
in the present paper this work is extended to 3d. The confining correlator 
D H (x) is expressed through the two-gluon gluelump Green's function, and 
for the latter one has to use the relativistic 3-body dynamics, which is well 
developed in Refs. [T3]-[TB] . Another part of our study is the calculation of 
the gluonic bound states: glueballs, gluerings and gluestars. We are using the 
same Hamiltonian approach as in the Ad case [T7| , [TS] to compute masses of 
these states in terms of the only parameter of the theory, a s (t). The resulting 
ratios M n / ^fa~ s are in good agreement with existing lattice data. 

The plan of the paper is as follows. In section 2 the main formalism is 
introduced, and the 3 body gluelump Green's function is calculated, together 
with a s (T). In section 3 the glueballs are calculated in the 3d gluodynamics, 
and in section 4 gluerings and gluestars are introduced and the masses are 
computed. In the last section a comparison to lattice and other approaches 
is given and a summary and prospectives are outlined. 



2 The large T limit of QCD 

The QCD action S at finite temperature 



S = S g + S q , S g = —J o dx 4 Jd 2 x(Fp 2 , S q = jf dx 4 Jd 2 x^(m-D)^ 

90 ° (1) 

at large T goes over into the effective 3d theory jHj-jOj, where the leading 
part of S g transforms as 

S 9 = ^fJ i F ik?d"x = i / (F*)Wx (2) 
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with g\ = g 2 T is a dimensionful coupling constant of the 3d theory. It is also 
known that quark degrees of freedom decouple while A 4 can be effectively 
integrated out at large T 0-0- 

Our aim is the calculation of field correlators in 3d theory using the 
gluelump Green's function, and to this end we absorb g 3 in F ik as Fi k = — F ik 
and we define the correlator in the same way as in the 4d case |l(J| ITT] 

— (tr(g 3 F ik (x)$(x,y)g 3 F em (y)®(y,x)) = D ikM x,y) = 

= (SuS km - 5 im 5 k i)D H (x - y) + ~[dihi5 km + perm}D±(x - y). (3) 

The spacial string tension a s is expressed through D H (x) in the usual 
way 013 EH 

° s = \\ DH ^ d2x - ( 4 ) 

Note that we keep index H in D H (x), appropriate for 4c? theory, since D H 
coincides with the high T limit of 4d magnetic correlator DjjL 4 . 

To proceed one can use as in [13] the Background Perturbation Theory 
(BPTh) [HI and write in the lowest order of BPTh (cf. Eq. (35) of [13]) 

D ",im( x iy) = -^H( tr a([ai,akMx,y)[ai,a m })) (5) 

c 

where $ and tr a are the parallel transporter and trace in adjoint representa- 
tion respectively. 

Writing [a«, a k ] = ia^a b k f abc T c one can define the 2-gluon gluelump Green's 
function G^ 29l \x,y) as 

D? kM {x,y) = ^tr a (r fe 7 de/ <(x)4(x)T c $(x, Z/ )T/af(y)a^( 2/ )) = 

g ^ N \~ l) {5ukm - 5 im 5 kl )G^ l \x, y) (6) 
and hence from (j3J) one has 

D H (x,y) = 9l{N l- 1) G^(x,y). (7) 
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The two-gluon gluelump Green's function G^ 2 ^ in 4c? was studied and 
numerically evaluated in ^3], and in what follows we shall repeat the same 
calculations for the lowest state in 3d. 

In the lowest order in g% the BPTh gives for G 2g ^ the following Fock- 
Feynman-Schwinger (FFS) path-integral representation [2~T] 

G^ l \x,y) = J o d Sl j o ds 2 (Dz^) xy (Dz^) xy e- K ^(W F ) (8) 

where 

W F = PP F exp ig 3 / Aidz\ 'expig 3 / A k dz\ >U F (x,y) (9) 

and 

U F (x,y) = exp 2^3 / dr 1 F(z l {T l )) e^>2ig 3 dr 2 F(z 2 (T 2 )) (10) 

JO JO 

As a next step we neglect the spin interaction of gluons due to the term 
[/p in (j§J), which was studied in ^3] and shown to be small (spin splittings 
< 10%), and consequently put U F — 1. The resulting Wilson loop in (JHJ), 
H^i?$(a;, y), consists of two gluon trajectories in W F and one straight-line 
trajectory in $(x, y). The vacuum average of this Wilson loop produces the 
area law factor which can be written as 

(W F $(x, y)) = exp(-a s (T)A min (x, y)) (11) 

where A min (x,y) = J c^xVdet g is the minimal area surface connecting the 
three trajectories discussed above. 

For small vacuum correlation length A, A hadron size, one can reduce 
the problem to that of a local Hamiltonian. |22[ |2"3] 

G^ l \x,y) = (x\e- H{23lh \y) (12) 

where H^ 2gl ^ is constructed in the same way, as for the qq case in [22]; and 
in the Ad case it was written in ^3] . One can deduce from (J12)) the equation 
for G^ 29 ^ (x , y) , separating 2d space coordinates and Euclidean time t 

-^-H^iv^G^iv^ ^Oj^Wn^^-rD (13) 
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In our case of d = 2 + 1, the corresponding 2d Hamiltonian has the form 
(see Appendix 1 of ^3] for details of derivation). 

= ^1 + ^ + ^ + *.{T){\n\ + |r 2 | + In - n|} (14) 

where /ij are einbein masses to be found from the minimization of eigenvalues 
22J; Pj, rj are 2d vectors of momentum and coordinate of the 2-th gluon, while 
the position of the adjoint source (the $>(x,y) trajectory) is chosen at the 
origin. 

To solve (fPfj) for the eigenvalues M^ 9l '(pii, H2) it is convenient to use the 
hyperspherical harmonics ^HJ EE! , with the global radius p, 

2 +2 , 2 n — r 2 n + n / 1C \ 
p = * +1, ^ = -JT^ i = -JT ] (15) 

and the Schroedinger-type equation for the lowest harmonics looks like (one 
can take p,\ = p.2 — P, as follows from minimization anyhow) 

_i_A 4 ^ + y(p)V = ^ (16) 

with 

d 2 3 dty r 
A4 = ^ + Jd~p' V{P) = J d0 * V ( Pl ' T & W 

One obtains 

V{p) = c v p, d, » ( 1 + -g- J (7 S = ca s (18) 

A good accuracy (around 1%) one obtains when using the simple variational 
wave function in (JTfij) 

vfr(p) = C exp (Vy ) , S(ai, v) = (5) + V l. (19) 



Minimization with respect to the variational parameter v yields 

3 (2 
2 U 



^ r i 5/2) l , ^, v*) = l( ~) (^r/5/2)) 2 / 3 , (20) 
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C = V2v 2 ,\y(0)\ 2 = C 2 = 2v*. 
Finally, the gluelump mass is obtained from minimization over // , 
+ E(fi, f ) = and nas 

m = M (no, u Q ) = 4ii = 4.52^. (21) 

The spectral decomposition of the gluelump Green's function has the usual 
form 

G< 2 ^(O,t;O,0) = G^ l \t) = £ ^fie-^l" (22) 

n=0 (2/^n) 

where |\l/ n (0)| 2 = §f = |fL We have used in the standard boson 
normalization ^7= = ^7= for each gluon, since /1 = /io is the average energy 

of the gluon, as one can deduce from (jl4j) . note that as in 
In this way one has from (j3J) 

ffs = = aw-P f (23) 

2 2 J 87r 

where £ is 

/oo /^y2 00 

= £ = £ £„. (24) 

One can see from (J2(JI21|) that the term with n = contributes to £ very 
little 

00 

£ = 0.063 + £ £ n . (25) 

n=l 

The sum over n in (J24|) is actually diverging, since G^ 2,9 ^ (t) as the Green's 
function of the Hamiltonian (|14|). should behave at small i as the free Green's 
function, i.e. G^ 2gl \t) ~ t~ 2 , and the integral / d 2 xG^ 29l \x) is logarithmically 
divergent at x = 0. This purely perturbative contribution to D H , is exactly 
cancelled by other correlators in the cumulant sum, as shown in |24j . One 
can show, in the same way as it was done in for the 4d case, that G^ 2gl ^ (t) 
has nonsingular behaviour after subtraction of the perturbative singularity 
0(t~ 2 ). In this derivation one should take into account that for small loops 
the area law behaviour (JTTJl is replaced by the quadratic area low, (W) ~ 

exp I ^ A min ) , see second ref. in [TU] for details. To estimate £, we 
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replace the linear interaction in H^ 2gV> by the oscillator one, (with lo ~ m /2 
to reproduce the asymptotics), and obtain 

r°° zdz 2cj . . 

^1 = (26) 

One can see, that £ is of the order of 1, but £ < 1 for reasonable values of Zq. 
Here introduction of zq regularizes the perturbative behaviour of G^ 2,9 ^ (x) at 
small x and one can see that this interplay of perturbative and nonpertur- 
bative regime at small x of vital importance for the resulting value of £, and 
hence of a s . This should be compared with lattice data, where one obtains 
for N c = 3 U 

v /^ = ^0.566 = 0.566Tg 2 (T); £ w « 1.00 (27) 



3 Glueball excitations in the 3d gluodynamics 

As was discussed above, the high-temperature QCD has the 3d gluodynamics 
as a limiting theory, and therefore it is interesting not only to compute field 
correlators and a s (T), as was done in the previous section, but also to find 
the spectrum of this limiting theory. 

In this section we shall calculate the spectrum of two-gluon (2g) and 
three-gluon {3g) glueballs in 3d SU(N) gluodynamics and compare results 
with lattice calculations. The lowest two-gluon glueball state was previously 
calculated in the framework of our formalism in ^2|- We define the creation 
operators for glueball states together with spin and parity J PC in Appendix 1. 
One can start from the general QCD string formalism with spin-dependent 
terms, as was done in ^7J ^] in the case of Ad QCD, but here we shall 
neglect both spin corrections and string rotation corrections [221 > important 
for L > 2, and consider the simplest version, calculating with accuracy of 
10% or better, so that one has the Hamiltonian for the 2g case 

P 2 

# = //+ — + a a r: HV n = E n {fj,)V n , (28) 

with the adjoint string tension a a = ^l^ si and the glueball mass is M n = 
En(fJ'o), where /io is the stationary point of E n {n), ^-\^ = ^ = 0. It is known 
[25, 26 1, that the accuracy of this procedure is around 5-7%. Instead of 
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computing (J28)) numerically, solving the 2d differential equation, we shall 
approximate solutions replacing a a r by the term (^7- + § ) and ^ is later 
to be found as well as [i from the stationary point equation. We shall call 
the resulting energy E n (fi, v) and corresponding masses M n . We expect (and 
calculations support this) that M n differ from M n within 10%. As a result 
one obtains a simple expression for M n 

M„ = 4^^, 71 = 1,2,3,... (29) 

Here n — L + 1, and in our approximation radial excited state + * is degen- 
erate in mass with the 2 + state. 

As was mentioned above, we disregard spin-dependent interaction and 
therefore the L = state corresponds to degenerate ++ ,2 ++ states while 
L = 2 state - to degenerate ++ , 2 ++ , 4 ++ states and so on. 

The connection ()29j) is universal for any SU(N) group, however to connect 
masses to fundamental spacial string tension a s one needs to use the relation 

2N 2 9 8 

° a = Js[2~[ asi ^ = ^ S (SU(3)), a a = -<7,(SU(2)) a a = 2^(SU(oo)) 

(30) 

The resulting ratios ^= are shown in Table 1 in comparison with lattice 
data from [2T] 1 . 

Table 1: The mass ratios M n j ^fo~ s obtained from the Hamiltonian 
vs lattice data for angular momentum L = 0, 1,2,3,4. For comparison two 
lowest masses obtained on the lattice [23 are given in parentheses. 

lr The spin-parity assignment in |27j is in general different from ours in Appendix 1, and 
therefore only states with reliably the same quantum numbers are shown in the Tables 1 
and 2. 
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IPC 
J 


n++ 9++ 
u , z 


n — v 9 — i- 

U , Z 


n++ 9++ 
u , z 

4 ++ 


n — i- 9 — V 
u , z 

4 -+ 


n++ 9++ 
u , z 

4++ ;6 ++ 


L 





1 


2 


3 


4 


Ml frT 

SU(2) 


4 fi1 
(4.72) 




7 Q8 
(7.82) 






SU(3) 


4.24 
(4.33) 


6 


7.34 
(7.15) 


8.48 


9.48 


SU(oo) 


4 

(4.065) 


5.65 


6.92 
(6.88) 


8 


8.94 



The 3g glueballs can be considered in the same way, as it was done in 
[T71 ITS] for the 4c? case. The Hamiltonian is similar to the one used in section 
2, Eq. (jTlj) and reads (we put /ii = /x 2 = /i 3 = /i). 

^ ) = ^ + ^^ + ^El^-r i |. (31) 

The hyperspherical component of the grand orbital momentum K satisfies 
an equation (see Appendix 2 for details of derivation) 



2/i I dp 



+ 



2K + 3 d 



P dp 



(32) 



with V(p) = cp, c = 2V2. 

As is well known (see e.g. ^H]) the solution of ()31|) is given with 1-2% 
accuracy by the equation 



(33) 



where 



W(p) 



2pp* 
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+ V(p), 



and p is the stationary point of W(p), W\p)\ p=po = 0. The mass M K , n = 
\ ^ + E Kn (p) \ is obtained to be 

M K , n = ^2^(3A) V4 ^ + ?n+^ 3/4 ; A ^ + ^ ^ + ^ (34) 

The resulting ratio of Mx,n/ are given in the Table 3 for the radial exci- 
tation number n = and K = L = 0, 1, 2, 3. 

Table 2: The three-gluon glueballs in the 3<i SU(N) gluodynamics. Mass 
ratios M^^j ^fo~ s are given for L = 0, 1, 2, 3. The lowest C = — 1 mass from 
lattice data [21] is given for comparison in parentheses (see footnote on the 
previous page). 



K=L 





1 


2 


3 


jPC 


1—3— 


1+-.3+- 


3 — , 1 — .5 — 


1+- 3+~ 


Ml,o 


6.01 


7.48 


8.68 


9.74 






(5.91,0—) 




(8.73,2—) 





Using Eq. one can also calculate radial excitations, which are ac- 

tually collective excitations corresponding to the collective variable p. The 
resulting mass ratios for n — 1, 2 are 9.39 and 12.35 respectively. 

One can now consider the pomeron and odderon trajectories, using mass 
formulas (J29)) and (J34)) . The pomeron trajectory corresponds to J = L + 2, 
where "2" comes from the sum of "gluon spins" as in the d — 3+1 case. In 
the d = 2 + 1 case (as well as for d = 3 + 1) the creation operators are of the 
form tr(D L / 2 E t D L / 2 E k ). 

From (|29p one has 

Mj = 8a a (L + 1) = 8a a (J - 1) (35) 

and the pomeron intercept is exactly unity, while the Regge slope corresponds 
to the relativistic potential case, since string corrections |2*2ll2*3*] are not taken 
into account. 
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The odderon trajectory is obtained similarly identifying J = L + 3, and 
finding from (|34|) (with accuracy better than 1%) analytic expression 

M 2 (J) = a(J- A odd ); (36) 

where a = 2 7 / 2 v / 3a s = 19.6a s ; A odd = 2-^- = 1.14. 

One can see, that surprisingly the odderon intercept appears higher, than 
that of pomeron, in contrast to the situation found in the d = 3 + 1 QCD 

EH- 



4 Gluerings and gluestars 

For number of gluons N > 3 one has two options in constructing gauge- 
invariant gluon configurations: 

a) gluon ring, consisting of N gluons, "sitting" on the closed loop of 
fundamental string, with the potential energy V r i ng = 07 J2i=j+i \ r i ~ r j\] 

b) gluon star, consisting of iV gluons connected by N adjoint string pieces 
with the string junction, where color indices are assembled in d abc (for N c = 3) 
or in higher order operators for N c > 3. The corresponding potential energy 
is V st ar = o"a Y^iLi \ r i ~ Y|, Y - is the string junction position. 

The masses of gluon rings and gluon stars can be computed in the hyper- 
spherical formalism using formulas given in Appendix 2. One has from 
(IATT3|) 

M { kn 9) , k 2 mm/ 2(n + l/2)\ 3/4 , s 

= 2(NA KN c 2 N )V 4 1 + y /wr - > , (37) 

with A kn =(k + N+%)(k + N-I), c jy = ^^l r(|), and 
f r f ng) = J^- • ^ ne can i mme diately see that gluestars are 14% 

M KN V CTs 

(N = 3) to 26% (N = 00) heavier than gluon rings and both saturate at 
large N, i.e. 2.24(A^ -> 00). 



A similar situation and the same ratios are obtained for D = 3 + 1 in 



5 Summary and discussion 

An important feature of our results is that our calculation does not use any 
parameters. Several approximations have been made in obtaining (J29j) . (J34)) : 
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1) the use of a local Hamiltonian, implying small vacuum correlation length 
A, A hadron size, and small mixing with hybrids due to the large mass 
gap [28J; 2) spin forces were considered weak and self-energy corrections van- 
ishing due to renormalization and gauge invariance (see appendix 4 of ^3] 
for discussion); 3) color Coulomb interaction between gluons was neglected 
due to the expected strong cancellation with the one-loop correlation, ob- 
served in Ad BFKL formalism (cf. discussion in [T71 ITS] ). In d = 4 glueball 
calculations [T7J [TH] this suppression of Coulomb resulted in good agreement 
of calculated and lattice measured glueball spectrum. 

Hence all our approximations are expected to be accurate within fifteen 
percent or better, which is supported by a good agreement with lattice data 
[77] in Tables 1,2. 

On theoretical side the d — 2+1 gluodynamics was considered also in |3TJ] . 
where the authors developed a Hamiltonian formalism for the calculation of 
a s and vacuum excitations. Their result for cr s corresponds to £ = 1 in our 
Eq. (|24|) and agrees with lattice data, however the glueball masses to the 
knowledge of the author are not yet available from the formalism of |3TJ] . 

Summarizing, we have obtained the confining field correlator, a s , gluelump 
and glueball masses for d = 2 + 1 SU (N c ) gluodynamics. Glueball masses for 
configurations of different types: simple glueballs for N = 2, and gluerings 
and gluestars for iV > 3 are computed explicitly without any parameters 
through a s . 

Glueball spectrum has a clear hierarchy in growing angular momentum 
L and in number of gluons N. In this and other respects the d — 2 + 1 
spectrum is similar to that of d — 3 + 1. 

The author is grateful to N.O. Agasian and A.B.Kaidalov for discussions. 
The work is supported by the Federal Program of the Russian Ministry of 
Industry, Science, and Technology No. 40. 052. 1.1. 1112, and by the grant for 
scientific schools NS-1774. 2003. 2. 
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Appendix 1 



Creation operators of the d = 2 + 1 SU(N) glueballs 

In d = 2 + 1 one can choose x 3 as the Euclidean time evolution coor- 
dinate, and denote Ei = F^, H 3 = F12. The gauge invariant creation 
operators are constructed from powers of E i) H 3) and higher operators also 
include powers of Di = di — igAi, with i = 1,2. The 2d angular momen- 
tum L and parity for a given creation operator is defined by the assign- 
ment of L — 1 to H 3 and D, L = to The 2d -parity is defined by 
PEi = —Ei, PH% = Hz,PDi = —Di. The combination DiEi vanishes 
due equations of motion. The C parity is defined from the transformation 
Ff k — > —(Fi k ) T , Df — > — -Df . The resulting operators together with L and 
J PC assignments are given in Table 3 for 2g glueballs and in Table 4 for 3g 
glueballs. 

Table 3: Creation operators for 2g glueball states Ei = F i3 , H 3 = F 12 in 
d = 2 + 1 SU(N) theory. 



L 




Operator 





0++ 







2++ 


tr(EiEk) symm. 


1 


2-+0-+ 


tr(EiH 3 ) 


2 


2++,0++ 


tr(H 3 H 3 ) 


2 


0++ 


tr(DiE k DiE k ) 


2 


2++ 


tr(DiE[D k Ei) symm , 


2 


4++ 


tr{P)iE k D[E rn ) symm. 


3 


2-+,0-+ 


tr[D t E k DiH 3 ) 


3 


4-+, 2-+ 


tr{DiE k DiH 3 ) symm. 


4 


2++0++ 


triDiH.DiH,) 


4 


6++, 4++, 2++ 


tr(DiH 3 D k H 3 ) sy mm. 


4 


6++ 


tr(DiD k EiD n D m E p ) symm. 
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Table 4: The same as in Table 4, but for 3g glueballs. 



L 


J PC 


Operator 




1 
1 
2 
2 
3 


1— 

3— " 
1+- 

3+-.1+- 
3— ,1— 
2-+, 4-+ 
l+-,3+- 


tr{EiE k E k ) 

t T ( £< j _£/j. .£// J symm. 

tr(E?H 3 ) 

tr((EiE k )H 3 ) symm. 

tr{EiH 3 H 3 ) 

triEiEkDtHs) 

tr{H 3 H 3 H 3 ) 



Appendix 2 



Calculation of the N— gluon glueball masses for TV > 3 

The Hamiltonian of iV-gluon glueballs has the form similar to (|31|). we 
fix the gluon einbein masses to be the same: /ii = /i 2 = ... = /i). 

MfJ, A 2 jV- 2 , \ / i n i\ 

g = ^-+ 2 ^ +y J v(ri,-r JV ) (A2.1) 

where Vjv can be of two different forms, for gluon rings (closed gluon chains 
with fundamental strings between adjacent gluons) 

N-l 

V { n ] = a s J2 l r m " r N = r (A2.2) 

i=0 

and for gluon stars ( a generalization of the F-type glueball for N c — 3 with 
the color symmetric matrix d abc at the string junction vertex, see ^H] for 
details) 

N 

V^^Elr.-Yl (A2.3) 
t=i 
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where Y is the position of the string junction, which for gluons can be taken 
to coincide with the cm. coordinate R = -k J2f=i r i- 

We approximate the eigenfunction of H (jA2.1|) for a given angular mo- 
mentum L by the hyperspherical component [?] 

<Mn, r 2 , ...r N ) = ** + % -Pk(vi, r 2 , ...t n ) (A2.4) 
P 2 

where p 2 = J2iLi( v i ~ an d Vk is the harmonic polynomial of Jacobi 
coordinates £ i; i = 1,N — 1, constructed of r 1; ...r N , 



v%+i) 



(A2.5) 



Averaging of over angular variables Q2N-2 in (2Af — 2) space 

of {^j} yields an equation for Xk(p) 

("27$ + Wkn ^ + If) MM = M Mp)xk(p) (A2.6) 
where Wkn{p) is 

^v(p) = ^ + 4 fc) (p), A^=(if + JV-|)(lf + J\r-|) (A2.7) 
and 

V 2 f)(p) = (y j f ) (r 1 ,r 2 ,...))n, /c = r, star. (A2.8) 

To calculate vjf (p) it is enough to know the angular average of one 
Jacobi vector, 



\( k )T(N - 1)T (-) 

V^ip) = a s NX^(\^\) u = a sP N X- = C&V. (A2.9) 



and A'-' -' = for k = r(ring), and \( star ^ = y^r-^ 2 -, as one can deduce 
from |AT5|) for i = N - 1. 

Finally the solution for Mku{p) can be found with ~ 1% accuracy from 
the minimum of Wkn(p) 

function of p jT^l UHl EH] • 



^(li) = W KN (p )+u,(n + ±), u 2 = -W"(p ) (A2.10) 
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1/3 

where p = ( ^ffeT ) > so that one has 



rr^Nip)__VN _ 3.1/3/^x2/3 A | 2(w+ 1/2) \ 



From ^M| M=W=0 one has 



^r"~ + ^' ^"la/vJ (A2 - 12) 



and finally the mass of the N gluon state is 



M KN _ on i/ 4 /. 2 ( n + 



3/4 



2£#* 1 + v /oA z/ (A2.13) 



where 

= NAkn(C^) 2 . 

For TV = 3 and K = 0, 5 3 = 18 and = 6.01,9.39,12.35 for n = 

0, 1, 2 respectively. 

For iV large, one can consider a gluon ring of A" gluons connected by 
fundamental string for any N c . 

Then asymptotically C KtN (N -> oo) « ^=j, Sjvioc = iV (f ) 1/4 and 



For iV = 3 this asymptotic formula yields ^!== ~ 6.72 which is not far 
from our exact value of 6.01. 

The gluon stars are obtained with A (star ) in cffi instead of \( rm 9\ which 

/ \(star) / / /V 1 \ ^"/^ 

yields for the mass the relative factor y fj^j = J^f- ■ i^r) ■ F° r N = 3 
this factor is 1.14, which means, that the 3-gluon star is 14% heavier than 
the 3-gluon ring. 
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